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1  Introduction,  Analysis,  and 
Application 


Background 

In  a  simple  towing-tank  experiment,  the  free  surface  excited  by  the  forward 
motion  of  a  barge  model  was  studied.  Sketches  of  the  model  and  tank  are  shown 
in  Figures  1-3.  Several  rather  important  observations  were  made  during  this  two- 
dimensional  (approximately,  since  the  model  spanned  the  channel)  study  of  a 
rectangular  barge  in  a  narrow  channel.  First,  when  the  towing  speed  was  low 
enough,  no  dispersive  surface  waves  were  created.  Instead,  the  free-surface 
elevation  at  the  bow  was  found  to  increase  as  the  model  increased  speed;  while, 
at  the  stem,  the  free-surface  drawdown  increased  with  speed.  This  latter  obser¬ 
vation  was  somewhat  obscured  by  wake  effects.  As  the  speed  of  the  body 
increased,  dispersive  surface  waves  were  created,  producing  rather  interesting 
reactions  on  the  model  in  the  narrow  channel.  These  reactions  will  be  addressed 
in  Phase  2  of  the  present  study.  In  this  report,  an  analytical,  three-dimensional 
added-mass  study  is  presented.  The  analytical  model  developed  herein  is  based 
on  the  low-speed  (0  to  4  fps  (0  to  1.2  m/s))  free-surface  behavior  observed  in  the 
tank  study. 


Analysis 

The  mathematical  model  developed  in  this  study  describes  a  barge  approach¬ 
ing  a  quay  wall  at  angle  a,  as  diagrammed  in  Figure  4.  In  that  figure,  we  see  an 
areal  water-plane  diagram  of  the  approaching  barge  on  the  left-hand  side  and  its 
mirror  image  on  the  right-hand  side.  The  use  of  the  mirror  image  allows  for  the 
effects  of  the  presence  of  the  quay  wall.  Two  inertial  coordinate  systems  are 
used,  along  with  three  linear  coordinates  {Xa,  Xc,  and  5)  on  the  physical  barge  and 
three  linear  coordinates  (Xa\  Xc\  and  s')  on  the  image  barge.  The  first  inertial 
system  (£,  £,  z)  has  c,  on  the  line  of  action  of  the  approaching  physical  barge;  the 
second  system  z)  has  L,'  on  the  line  of  action  of  the  barge  image.  The 

origins  of  the  two  inertial  coordinate  systems  are  at  the  same  point  on  the  quay 
wall,  where  the  calm-water  level  intersects  the  wall.  The  relationships  among 
these  are 

£'  =  £cos(2a)  +  ^sin(2G!)  (la) 
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Figure  2.  Plan  view  sketch  of  model  and  tank  (not  to  scale) 
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Figure  3.  Cross-sectional  view  sketch  of 
model  and  tank  (not  to  scale)  (to 
convert  inches  to  millimeters, 
multiply  by  25.5) 
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Figure  4.  Water-plane  sketch  of  barge  approaching  the  quay  and  its  mirror  image 


sin  (2a)  +  {  cos  (2a) 

(lb) 

Referring  to  the  sketch  in  Figure  5,  the  relationships  between  the  physical  barge 
coordinates  are 

(2a) 

x,=s-j+x0=f  +  m  +  x0 

(2b) 

Xc  =J~— +X0  =£+  L  +  C(t)  +  X 0 

(2c) 

where 

L  =  wetted  length  of  the  barge 

i  ( t )  =  distance  between  the  leading  edge  of  the  barge  along  the  line  of 
action  in  the  center  plane  of  the  barge  and  the  quay  wall 

Xo  =  distance  between  the  nodes  of  the  tree-surface  deflection  and 
intersections  of  the  bow  and  stem  with  the  calm-water  plane  along 
the  center  plane 
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Figure  5.  Elevation  sketch  of  physical  barge  and  free  surface 

Between  the  barge  coordinates  of  the  physical  and  image  barges,  the  following 
relationships  exist: 

W 

xa,c  =  xa,c  [2cos  (a)  - 1]  + — sin  (2a)  +  2Tacsin  (a) 

Ca,c  =  Ca,c [2cos(ar)  - 1  ]  -  Xa  c  sin  (2a)  +  Wcos*  (a) 

+  2Yoccos(a) 

where  the  double  subscripts  (a,c)  are  used  for  brevity.  Referring  to  Figure  4,  the 
distances  YJj)  and  Y^t)  [=  Ya(t)  +  L  sin  (a)]  are  those  between  the  respective  bow 
and  stem  barge  comers  on  the  starboard  side  and  the  wall  on  the  water  plane.  We 
should  note,  here,  that  in  the  application  of  the  analysis,  two  barge  orientations  are 
of  special  interest.  The  first  is  when  the  barge  is  traveling  parallel  to  the  quay  wall  at 
a  separation  distance  Ya  from  the  wall  to  the  barge  side.  For  this  condition,  a  =  0, 
and  the  image  barge  coordinates  are 

K=XaP=K  (2f) 

Cap  =Cap+W  +  2Yacp=C'p  (2g) 


(2d) 

(2e) 


where  W  is  the  barge  wetted  width  and  the  subscript  p  refers  to  the  parallel  travel. 
The  second  orientation  is  when  the  barge  is  normal  to  the  wall,  and  for  this,  the 
subscript  n  is  used.  The  image  barge  coordinates  for  this  orientation  are 


K=-xm+2  n 

c  =-c=-c 


(2h) 

(2i) 


Chapter  1  Introduction,  Analysis,  and  Application 


For  the  general  case  considered  here,  we  can  assume  that  the  barge  has  a 
constant  acceleration  (or  deceleration),  or  that  the  barge  is  running  at  a  constant 
speed.  Hence,  the  barge  speed  is 

U  (t)  =  U0+  a0t  (3) 

where  Uo  is  an  initial  barge  speed  and  ao  is  a  constant  acceleration  if  positive,  or  a 
constant  deceleration  if  negative.  The  origins  of  Xa  and  Xc  are  at  the  free-surface 
nodes  of  the  respective  bow  and  stem  wave  forms,  as  sketched  in  Figure  5. 


Hydrodynamics  about  the  physical  barge 

The  analysis  is  first  applied  to  the  physical  barge  (as  opposed  to  the  image 
barge).  Our  attention  is  focused  on  the  free  surface  in  the  neighborhoods  of  the 
bow  and  the  stem.  From  experimental  observations,  the  free  surface  is  assumed 
to  have  bow  and  stem  profiles,  as  sketched  in  Figure  5.  Ahead  of  the  bow,  there 
is  a  crest  at  s  =  L/2  +  6  and  a  node  at  5  =  L/2  -  X0.  At  the  stem,  a  trough  is  at  5 
=  -L/2  +  5,  and  the  node  is  at  s  =  -L/2  -  X0.  The  distance  8  is  to  be  determined 
by  rake  angles  at  the  bow  and  the  stem.  In  both  -Xq  <Xa<Xo  and  -Xo  Xc  <  Xo, 
the  free  surface  is  assumed  to  have  the  properties  of  a  ship-generated,  deep¬ 
water,  linear  wave. 

From  McCormick  (1973),  the  “half-wavelength”  in  these  regions  is 

2  X0=-U2  (4) 

g 

To  satisfy  the  assumed  free-surface  profile,  the  following  potential  function  is 
assumed: 


<PpB  ~  Vbow  +  stem 


(5) 


(PpB  -  “^0) 


1 


1 


1 


cosh  (AX  )  cosh  (BO  cosh  (AXC)  cosh  (BO 


cosh  [C(z+h)\ 


cosh  (Ch) 


where  the  coordinate  is  fixed  to  the  barge  and  the  subscript  PB  refers  to  the 
physical  barge.  The  £-term  results  from  the  assumption  that  there  is  a  negligible 
beam- wise  velocity  component  over  the  beam  of  the  barge  traveling  slowly.  That 
velocity  component  has  nonzero  values  beyond  the  barge  sides.  The  wave  form  is 
traveling  with  respect  to  the  inertial  coordinate  system,  since  £  +X0  in  Equa¬ 
tions  2  and  3  is  a  function  of  time.  The  coefficient  0(0  varies  in  time  due  to  the 
time-dependent  velocity  of  the  barge.  Because  the  velocity  is  assumed  to  be 
small,  the  coefficient  is  a  weak  function  of  t.  By  this,  we  mean  that  the  time 
derivatives  of  the  coefficient  are  of  second  order.  The  first  terms  in  the  brackets 
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of  Equation  5  represent  the  bow  wave-form  potential,  while  the  second  terms 
represent  that  of  the  stem. 

We  assume  that  the  heights  of  the  bow  and  stem  wave  forms  are  equal.  Note  that 
when  viscosity  is  considered,  the  stem  wave  form  is  reduced  in  height  due  to  the 
vortices  being  shed  from  the  bottom.  This  is  discussed  later  in  the  report. 

To  obtain  the  shape  of  the  free  surface  associated  with  the  velocity  potential,  we 
apply  the  linearized  Bernoulli’s  equation  at  the  free  surface.  The  resulting 
expression  for  the  free-surface  displacement  is 

tanh  (AXa )  1  tmh(AXc)' 

cosh  (AXa  )  cosh(Z?£")  cosh  (AXC ) 


cosh  (B£)  J 

Figures  6  and  7  show  plots  of  the  free-surface  displacement  described  by 
Equation  6.  In  the  region  of  the  barge,  the  free  surface  is  replaced  in  the  plot  by 
the  boundary  surface  of  the  barge.  The  free-surface  form  approximately  predicted 
by  expression  in  Equation  6  is  supported  by  experimental  observations.  Note  that 
for  a  time-dependent  speed,  the  separation  length  between  the  origin  of  Xa  and 
the  wall  can  be  represented  by 

(t  +  XQ)  =  -\udt  (7) 


---3 


AOU 


r=0 


Hence,  we  can  also  write  d(  £  +  X0)/dt  =  -U(t).  The  approximation  in 
Equation  6  is,  again,  based  on  the  assumption  that  the  barge  moves  slowly 
enough  to  cause  the  coefficients  <1>,  A,  and  B  to  be  weak  functions  of  time.  By 
this  we  mean  that  the  time  derivatives  of  these  functions  are  of  second  order. 

In  the  vertical  direction,  the  behavior  of  the  velocity  potential  of  Equation  5 
must  be  like  a  linear  wave.  This  behavior  satisfies  the  bottom  boundaiy  condi¬ 
tion,  which  requires  that  there  be  no  normal  flow  at  z  =  -h.  The  parameter  C  is  the 
wave  number  for  this  linear  wave.  The  wavelength,  X,  is  estimated  from  Equa¬ 
tion  4,  so  that 


2n  _  2n  _  g 


(8) 


At  the  wave  crest  at  the  bow  and  trough  at  the  stem  (Xa  =  X0  +  b  and  Xc 
=  X0  +  5,  respectively),  we  require  that  the  slope  of  the  free  surface  be  equal  to 
zero.  Applying  this  condition  along  the  center  plane  of  the  barge,  where  £  =  0, 
one  finds  that 
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Figure  6.  Free  surface  surrounding  jumbo  open-hopper  barge  (JOHB)  at  a 
speed  of  1 0  fps 


Figure  7.  Free  surface  surrounding  bow  of  JOHB  at  a  speed  of  10  fps 
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dlj 


PB„ 


■  + 


tanh2  ( AXa ) 


_  A*®U  -1 

X.=xt+s,c= 0  g  Lcosh3(^)  '  cosh (AXa ) 


=  0 


(9) 


The  nontrivial  solution  of  Equation  9  is 
sinh2  [A(X0  +  <?)]- 1  =  0 


(10) 


Using  Equation  4  to  substitute  for  X0,  substituting  for  6  using  the  geometry  of 
Figure  5,  and  solving  for  the  parameter^  yields 


A  = 


L7627g 


U2\7T  +  - 


1 


(11) 


tan# 


The  velocity  potentials  at  the  bow  and  stem  should  satisfy  the  equation  of 
continuity  in  the  form  of  Laplace’s  equation 

v>  =  0  (12) 

Applying  Laplace’s  equation  to  the  velocity  potential  of  Equation  5  yields 


V2<PPB  =  -^/)COSh[C(,Z  +  /7)]^2 


B 2 


cosh  (Ch) 

1 


+C2 


cosh(  AXa ) 
1 


1 


-1  |  tanh2 (AXa) 

cosh3(AYa)  cosh(AY0) 

tanh 2  (BO 


1 


cosh(fi^) 


cosh3  (5^)  cosh(S^) 
1 


(13) 


cosh(^Ya)  cosh(fi^) 


■  =  0 


On  the  center  line  of  the  barge  in  the  region  of  the  bow,  Equation  13 
simplifies  to 


PB„ 


1  =  A2 

-1  (  tanh2(yfX0) 

-  **■ 

k=o 

cosh3(AYa)  cosh(J4Ara) 

52 


-  +  ■ 


C2 


(14) 


cosh  (AXa )  cosh  (AXa ) 
Letting  Xa  go  to  zero  at  the  bow  yields 
-A2-B2  +C2  =  0 


=  0 


(15) 
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Solving  for  B, 


B  =  Jcr^A2 


(16) 


Next,  we  know  that  points  of  the  maximum  and  minimum  free-surface 
displacements  from  the  calm-water  plane  at  the  barge’s  center  plane  (£  =  0)  are 
the  points  of  the  water’s  maximum  energy  absorption  from  the  barge.  The 
potential  energy  per  unit  volume  of  the  water  at  these  points,  which  is 
proportional  to  the  free-surface  displacement  at  these  points  (rioo  and  t|co  in 
Figure  5),  is  equal  to  the  water’s  kinetic  energy  per  unit  volume  at  the  bow. 
Hence,  by  incorporating  the  expression  in  Equation  6,  we  can  write  the  following 
expression  for  the  maximum  free-surface  displacements: 

„  =_„  =  A$>U  [  tanh  [A(X0  +  S)]  \  ~  \  Jj2 
K  ,Co  g  \cosh[A(X0  +  S)]}  2  g 


where  O  is  to  be  determined.  Recall  thatX0  =  %U2llg  from  Equation  4.  Also, 
from  the  geometry  shown  in  Figure  5,  8  =  T|0(/tan  (0).  The  result  in  Equation  17 
can  be  rearranged  to  obtain  the  following  expression: 


<D  = 


U 


cosh2 


Air 

2  g 


7U  +  - 


1  \ 


tan# 


2  A 


sinh 


AU2(  1 

JC  +  - 


2  g 


V 


tan#. 


(18) 


where  0  is  the  rake  angle. 

In  summary,  the  respective  coefficients  A,  B,  C,  and  O  of  Equations  11,  16, 

8,  and  18  are  parametric.  The  parameters  are  time  dependent  if  the  motion  of  the 
barge  is  unsteady,  i.e.,  if  U  =  U(t). 

In  the  following  section,  the  velocity  potential  and  free  surface  of  the  image 
barge  is  presented.  This  is,  again,  done  in  the  barge  coordinate  system.  Hence,  we 
use  the  coordinate  relationships  in  Equations  2d  and  2e  in  expressions  of  the 
forms  of  Equations  5  and  6. 


Free  surface  about  the  mirror-image  barge 

The  velocity  potential  and  free-surface  displacement  expressions  for  the 
mirror-image  barge  are  similar  to  those  in  Equations  5  and  6,  respectively.  The 
differences  are  in  the  coordinates,  as  expressed  in  Equation  1.  In  terms  of  the 
physical  barge  coordinates,  the  velocity  potential  for  the  mirror  image  is 
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_ 1 _ 1 

cosh(/fZ^)  cosh(fif*) 


_ 1 _ 1 

cosh(^.)  coshf/if) 


cosh[C(z+  /»)] 
cosh(C7?) 


(19) 


and  the  free-surface  displacement  is 


V  IB 


A<bU 


g 


tanh  (AX'a ) 


iw\h(AX’c ) 


cosh(^Y')  cosh(Z?£")  cosh(AX'c )  cosh(B£") 


(20) 


where  Za  c  and  t,a.c  are  obtained  from  Equations  2d  and  2e,  respectively. 


The  results  in  Equations  5  and  20  can  be  added  together  to  obtain  the 
expression  for  the  velocity  potential  at  any  point  in  the  physical  flow  field  (£,  <  0 
in  Figure  4).  That  is,  for  any  point  in  the  physical  flow  field,  the  velocity  can  be 
obtained  from  the  following  velocity  potential: 


(p  =  (pPB  +  (pm  (21) 

In  addition,  the  free-surface  displacement  expressions  of  Equations  6  and  21 
are  added  together  to  obtain  the  free-surface  geometry  in  the  physical  field.  The 
result  is  that 


7l  =  r?PB+  Vm  (22) 

We  shall  use  these  expressions  to  obtain  the  added  mass  excited  by  the  physical 
barge. 


Kinetic  energy  and  added  mass 

As  derived  by  Karamcheti  (1966)  and  other  books  dealing  with  advanced 
fluid  mechanics,  the  kinetic  energy  of  a  potential  flow  that  is  caused  by  a  body 
traveling  at  a  speed  U is  mathematically  expressed  as 

=  (23; 

s  s 

where  (referring  to  the  sketches  in  Figures  8-10), 
p  =  mass  density  of  the  ambient  water 

n  =  normal  unit  surface  vector  at  any  place  on  the  surface  area  (S)  of 
the  barge 

mw  =  added  mass 
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Figure  8.  Relative  and  absolute  flow  directions  in  the  center  plane 
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The  added  mass  is  a  measure  of  the  inertial  reaction  of  the  fluid  on  the  barge. 
From  Equation  23,  then,  the  added-mass  expression  is 

=  (24) 

s  s 

Equation  24  is  solved  numerically.  The  results  are  applied  to  a  jumbo  open- 
hopper  barge  for  various  values  of  the  parameters  and  approach  angle. 

Appendix  A  presents  the  computer  program  used  to  obtain  the  results  described 
herein.  Example  input  and  output  runs  are  shown  in  Appendix  B. 


Application  to  the  Surfaces  of  the  Physical  Barge 

To  determine  the  added  mass,  one  applies  Equation  24  to  each  surface  of  the 
body  using  each  of  the  velocity  potentials.  On  the  bow  surface, 


dcp 

dn 


bow 


=  ^sin#  ~^r-cos6 


% 


dz 


(25) 


while  on  the  stem, 


d(p 

dn 


stem 


.  „  d(p 
sin# — r-cos# 
oz 


(26) 


as  seen  in  Figure  10. 
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Results 


Before  describing  the  results,  a  brief  discussion  should  be  given  regarding 
the  nature  of  the  added  mass.  As  previously  mentioned,  the  added  mass  is  a 
measure  of  the  inertial  reaction  of  a  fluid  on  the  motion  of  a  body.  When  the 
body  is  at  rest,  the  added  mass  is,  naturally,  equal  to  zero.  When  a  body  travels  at 
a  constant  velocity,  there  is  a  convective  acceleration  of  the  fluid.  Hence,  even  at 
a  constant  body  speed,  there  is  an  inertial  reaction  from  the  fluid,  as  exhibited  by 
the  added  mass.  When  the  body  accelerates,  the  total  inertial  reaction  equals  the 
virtual  mass,  which  is  the  body  mass  plus  the  added  mass.  The  force  on  the  body, 
then,  is  the  product  of  the  virtual  mass  and  the  body’s  acceleration. 

Now,  concerning  the  results  of  this  study,  Equation  24  is  applied  to  the 
jumbo  open-hopper  barge  (JOHB),  for  the  most  part.  A  comparison  of  one  of  the 
parametric  conditions  for  the  JOHB  is  made  with  the  standard  open-hopper  barge 
(SOHB).  Since  the  geometries  of  the  two  barges  are  similar,  the  behaviors  with 
the  parameters  are  also  similar.  For  the  two  barges,  the  values  listed  in  Table  1 
apply. 


Table  1 

Comparison  of  Parametric  Conditions,  Jumbo  and  Standard 
Open-Hopper  Barges 

Parameter 

JOHB 

SOHB 

Wetted  length  (L),  ft  [m] 

195  [59] 

175  [53] 

Wetted  beam  (W),  ft  [m] 

35  [11] 

26  [8] 

Draft  (D),  loaded,1  ft  [m] 

101  [3] 

101  [3] 

Draft  (D),  unloaded,  ft  [m] 

3  [0.9] 

3  [0.9] 

Displacement,  loaded,  short  tons  [million  kg] 

1,500  [1.4] 

1,000  [0.9] 

Displacement,  unloaded,  short  tons  [million  kg] 

279  [0.25] 

179  [0.16] 

Rake  angle,  e  (for  computations) 

Variable 

Variable 

1  The  loaded  draft  is  used  as  the  standard. 

Effect  of  Barge  Speed  in  Open  Waters 

Since  the  barge  speed  is  assumed  to  have  low  values  in  confined  waterways, 
the  maximum  speed  value  used  herein  is  4  fps  (1.2  m/s).  The  “added  weight” 


(m»S)  is  presented  in  Figure  1 1  as  a  function  of  barge  speed  for  the  JOHB 
traveling  in  open  waters,  where  the  operating  water  depth  (h)  is  18  ft  (5.5  m).  The 
JOHB  parametric  values  are  listed  in  Table  1;  a  rake  angle  of  45  deg  was  used  for 
the  JOHB  in  all  computations.  The  added  weight  is  chosen  so  that  a  direct 
comparison  can  be  made  with  the  displacement  tonnage  of  the  barge.  The  “jump” 
in  the  plot  could  result  from  the  passage  of  the  barge  through  the  critical  speed  in 
the  wave-resistance  curve  (see,  for  example,  McCormick  1973).  Further 
investigation  of  this  is  warranted  in  Phase  2  of  this  research. 


Effect  of  Water  Depth  in  Open  Waters 

In  Figure  12,  the  added  weight  is  shown  as  a  function  of  water  depth  at  a 
position  well  away  from  the  wall.  Again,  the  parameters  used  are  those  of  the 
JOHB  in  Table  1  (rake  angle,  45  deg). 


Effect  of  Barge-Wall  Separation  and  Barge  Load 

Because  of  the  rake  angle  of  the  bow  (45  deg),  the  only  separation  effect 
occurs  when  the  barge  travels  parallel  to  the  quay  wall.  The  wall  effect  extends 
only  several  feet  from  the  wall,  as  can  be  seen  in  the  results  presented  in  Figure 
13.  Hence,  for  all  but  the  smallest  of  approach  angles  (a),  there  is  no  wall  effect, 
according  to  the  theory.  This  is  due  to  the  low  barge  speed.  The  results  are 
presented  for  the  JOHB  under  both  loaded  and  unloaded  conditions. 


Rake  Angie  Effect 

Although  the  rake  angles  of  both  the  JOHB  and  the  SOHB  are  fixed,  the 
effect  of  rake  angle  is  presented  in  Figure  14.  As  can  be  seen  in  that  figure,  there 
appears  to  be  an  angle  for  which  the  added  weight  is  maximum.  The  purpose  of 
this  exercise  is  to  help  barge  designers  determine  an  optimum  configuration. 


Comparison  of  the  JOHB  and  SOHB 

In  this  case,  the  JOHB  and  SOHB  are  compared  under  loaded  conditions, 
operating  in  1 8  ft  (5.5  m)  of  water,  at  a  speed  of  3  fps  (0.9  m/s).  For  the  JOHB, 
the  added  weight  is  3.44  tons  (3,120  kg)  and,  for  the  SOHB,  2.56  tons  (2,322  kg). 
Although  the  scale  factors  for  the  two  barges  differ  for  both  length  and  width,  the 
behaviors  in  Figures  1 1-14  can  be  expected  to  be  the  same  for  the  SOHB. 
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Figure  1 1 .  Added-weight  variation  with  speed  of 
the  loaded  JOHB  traveling  parallel  to 
the  wall  in  18-ft-deep  open  waters 


Figure  12.  Water  depth  effect  on  the  added  weight 
of  the  loaded  JOHB  traveling  parallel  to 
the  quay  wall  at  a  speed  of  3  fps 
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Added  weight  (tons) 


— ■ —  Draft  =  10  ft 
-0- Draft  =  3  ft 


Ya  (ft) 


Figure  1 3.  Wall  effects  on  the  added  weight  of  the  loaded  and  unloaded  JOHB 
traveling  parallel  to  the  quay  wall  in  18  ft  of  water  at  a  speed  of 
3  fps.  Here,  Ya  is  the  distance  between  the  starboard  comer  of  the 
barge  and  the  wall,  measured  perpendicular  to  the  wall) 
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Figure  14.  Effect  of  rake  angle  of  the  JOHB  operating  in  18-ft-deep  open 
waters  at  a  speed  of  3  fps 


3  Discussion  and  Conclusions 


The  objective  of  this  study  is  to  mathematically  model  the  added  mass  of  a 
barge  tow  operating  in  restricted  waters,  i.e.,  waters  of  finite  depth  that  are 
bounded  by  a  vertical  wall.  The  model  is  based  on  the  free  surface  altered  by  a 
box-lighter  (a  barge  with  vertical  sides,  bow,  and  stem)  towed  in  the  wave/ 
towing  tank  at  the  Johns  Hopkins  University. 

The  barge  model  was  towed  using  a  gravity  towing  system.  Hence,  the 
acceleration  in  the  ideal  case  is  constant.  The  realities  of  viscosity  and  the 
presence  of  a  vertical  channel  wall  placed  in  front  of  the  towed  model,  however, 
produced  a  varying  resisting  force  and  therefore  a  varying  acceleration  over  the 
test  runs.  One  cause  of  the  varying  acceleration  was  the  changing  added  mass,  as 
is  discussed  later. 

The  free-surface  deflections  resulting  from  the  accelerating  barge  can  be  seen 
in  Figures  15-17,  where  the  respective  accelerations  vary  from  the  highest  to  the 
lowest  value.  Because  both  the  bow  and  stem  of  the  model  were  vertical,  the 
crest  of  the  bow  wave  and  the  trough  of  the  stem  wave  (due  to  drawdown)  were 
at  their  respective  barge  faces.  The  drawdown  is  most  apparent  in  Figure  17.  The 
flow  separation  at  the  transom  of  the  model  absorbed  much  of  the  energy  trans¬ 
mitted  to  the  water.  This  resulted  in  a  shallower  trough  than  would  be  observed 
under  ideal  (inviscid)  conditions.  The  free-surface  profile  at  and  near  the  bow 
could  be  seen  in  the  water  between  the  barge  and  the  tank  wall.  The  clearance 
between  the  two  structures  was  approximately  0.2  in.  (5  mm).  The  side 
“channels”  were  regions  of  relatively  high  viscous  losses. 

The  predicted  effect  of  the  barge  speed  on  added  weight  can  be  seen  in 
Figure  1 1  for  the  JOHB  operating  in  open  waters  where  the  water  depth  is  18  ft 
(5.5  m).  For  an  unexplained  reason,  the  added  weight  experiences  a  “jump”  in 
values  between  2  and  3  fps  (0.6  and  0.9  m/s).  For  speeds  higher  than  4  fps 
( 1 .2  m/s),  the  values  of  the  added  weight  drop  as  the  speed  increases.  This  is 
attributed  to  fact  that  the  theory  is  a  low-speed  theory. 

The  effects  of  water  depth  on  the  added  weight  are  shown  in  Figure  12  for 
the  JOHB  traveling  at  3  fps  (0.9  m/s)  in  18-ft  (5.5-m)-deep  open  waters.  For  the 
depths  studied,  the  increase  in  added  weight  with  increasing  depth  is  nearly 
rectilinear. 
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Figure  15.  Free-surface  profile  amidship  for  a  rectangular  model— high 
acceleration 


Figure  16.  Free-surface  bow  to  stern  for  a  rectangular  model— moderately  high 
acceleration 
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Figure  1 7.  Free-surface  profile— moderate  acceleration 


In  Figure  13,  the  effect  of  the  separation  distance  ( Ya )  between  the  quay  wall 
and  the  JOHB  traveling  parallel  to  the  wall  is  shown.  It  can  be  seen  in  that  figure 
that  the  wall  effect  at  a  3-fps  (0.9  m/s)  barge  speed  is  observed  only  within 
several  feet  from  the  wall,  for  both  loaded  and  unloaded  conditions. 

The  effect  of  varying  the  rake  angle  (0)  of  the  bow  and  stem  of  the  JOHB  is 
seen  in  Figure  14.  In  that  figure,  we  observe  a  peak  value  of  the  added  weight, 
occurring  between  angles  of  50  and  55  deg.  Depending  on  the  added-weight 
effect  on  performance,  this  curve  can  be  used  to  somewhat  optimize  the  system 
by  selecting  a  0-value  that  is  both  practical  and  corresponds  to  the  lowest  inertial 
resistance. 

Figure  1 8  shows  a  plot  of  the  velocity  potential  contours  in  the  water  sur¬ 
rounding  the  JOHB  at  a  velocity  of  3  fps  (0.9  m/s),  as  predicted  by  Equation  5.  A 
detail  of  the  flow  around  the  bow  of  the  barge  is  shown. 

Finally,  during  the  experiments  on  the  model  barge  motions  in  the  wave/ 
towing  tank,  the  barge  acceleration  was  measured.  The  acceleration  was 
integrated  twice  to  obtain  the  velocity  and  the  displacement  of  the  model.  In 
these  tests,  a  vertical  barrier  was  alternately  placed  in  the  path  of  the  model  to 
simulate  the  presence  of  a  wall.  The  results  of  these  tests  are  presented  in 
Figures  19-22.  Although  a  number  of  towing  weights  were  used,  only  two  are 
represented  in  those  figures.  For  the  low  towing  weight  (lower  acceleration  runs), 
the  presence  of  the  wall  had  little  effect  on  the  speed  and  acceleration,  as  can  be 
seen  in  Figure  19,  where  t  (/)  is  finite  (as  defined  in  Figures  5  and  8),  and  in 
Figure  20,  where  t  (/)  is  infinite.  For  the  high  towing  weight  (higher 
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Figure  18.  Velocity  potential  contours  for  JOHB  at  a  speed  of  3  fps 

acceleration  runs),  there  is  a  pronounced  decrease  in  the  acceleration  and  the 
corresponding  velocity.  Since  the  acceleration  is  a  measure  of  the  inertial 
reaction  of  the  water  on  the  model,  we  can  conclude  that  the  wall  causes  an 
increase  in  the  added  mass  for  the  higher  accelerations.  This  effect  is  also 
partially  produced  by  the  buildup  in  the  water  between  the  bow  and  wall  at  the 
higher  accelerations. 

For  accelerations  resulting  in  speeds  that  are  not  “low,”  another  phenomenon 
is  observed.  That  is,  a  traveling  wave  is  produced  by  the  initial  bow  motion  that 
travels  to  the  wall,  reflects,  and  returns  to  the  model,  bringing  the  model  to  a 
virtual  halt.  This  effect  is  actually  beneficial  for  barge  tows  that  are  traveling  at 
higher  speeds  toward  vertical  barriers,  since  it  helps  to  prevent  full-speed 
collisions  with  the  wall. 

This  preliminary  analysis  has  several  limitations,  and  further  study  is  recom¬ 
mended.  The  form  of  the  velocity  potential  of  Equation  5  is  based  on  experi¬ 
mental  observations.  Although  Figures  6  and  7  indicate  that  this  potential 
function  is  a  reasonable  prediction,  the  governing  equation  for  potential  flow 
(Laplace’s  equation)  is  met  only  at  specific  points  in  the  fluid  domain.  Further¬ 
more,  the  boundary  conditions  for  this  boundary-value  problem  are  met  only  at 
specific  points,  not  the  whole  boundary  surface  of  the  barge.  A 
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Figure  19.  Model  displacement,  velocity,  and  acceleration  under  lower  force  without  the  wall  in  place 


Figure  20.  Model  displacement,  velocity,  and  acceleration  under  lower  force  with  the  wall  in  place 


22 


Chapter  3  Discussion  and  Conclusions 
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Figure  22.  Model  displacement,  velocity,  and  acceleration  under  higher  force  with  the  wall  in  place 
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boundary-element  method  analysis  of  the  problem  is  planned,  and  will  solve  the 
complete  boundary- value  problem. 

In  conclusion,  the  analysis  presented  herein  is  primarily  heuristic,  being 
based  on  experimental  observations.  The  assumptions  made  limit  the  analysis, 
particularly  in  predicting  the  values  of  the  added  mass  or  added  weight;  however, 
it  is  believed  that  the  analysis  does  predict  the  behavior  of  the  added  mass  with 
varying  parameters.  A  follow-on  study  will  improve  the  quantitative  predictions, 
since  that  study  will  incorporate  the  boundary-element  method,  which  requires 
few  of  the  assumptions  made  herein. 

Finally,  the  analysis  presented  herein  is  considered  to  be  superior  to  the 
potential  theories  that  were  based  on  two-dimensional  models,  as  in  the  report  by 
Wendel  (1956)  and  the  report  by  Brennen  (1982). 
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Appendix  A 

FORTRAN  Computer  Code 


This  appendix  presents  the  computer  program  that  was  used  to  obtain  the 
results  described  in  this  report.  The  added  mass  of  the  barge  is  calculated  using 
Equation  24.  A  fourth-order  Runge-Kutta  technique  is  used  to  perform  the 
numerical  integrations  over  the  bow  and  stem  surfaces  of  the  barge. 

Note  that  the  program  uses  a  substitution  for  the  function  l/cosh(x),  which 
occurs  in  the  velocity  potential  function  and  its  derivatives.  Large  argument  in 
the  cosh(x)  function  cause  errors.  The  hyperbolic  secant  is  actually  desired,  but 
FORTRAN  does  not  include  this  function  intrinsically.  To  avoid  this  problem, 
the  trigonometric  identity 

tanh2  (x)  +  sech2  (x)  =  1  (Al) 

is  used  to  substitute  for  l/cosh(x),  so  that 


- =  sech(x)  -  -y/l  -tanh2  (x)  (A2) 

cosh(x) 

Since  tanh(x)  is  well  behaved  as  x  goes  to  infinity,  this  substitution  makes  the 
program  much  more  reliable. 

Use  of  the  program  is  self-explanatory.  Upon  execution,  it  prompts  the  user 
for  the  necessary  input  defining  the  barge  geometry  and  the  test  conditions. 
Modification  of  the  dimensions  of  the  barge  requires  modification  of  the 
computer  code.  The  appropriate  variables  are  annotated  in  the  code. 


PROGRAM  BargeAddMass 

C  This  program  is  used  for  the  WES  barge  added-mass  study.  There 
C  are  two  barges  of  interest.  The  first  is  the  Jumbo  Open-Hopper 
C  Barge,  195  feet  in  length,  and  the  Standard  Open-Hopper  Barge, 
C  175  feet  in  length.  Only  the  Jumbo  Barge  will  be  considered 

C  here.  The  integrations  for  the  added-mass  will  be  done  by 

C  a  fourth-order  Runge-Kutta  technique: 
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no  o  o  o  o  o 


INTEGRAL  =  deltaX/3  [f(Xl)  +4f(X2)  +2f(X3)  +4f(X4)  ... 

...  +2f(Xn-l)  +4f(Xn)  +f(Xn+l)  ] 

INTEGER  j,k 

DOUBLE  PRECISION  pi,  g,  rho,  degrees2rads 
DOUBLE  PRECISION  L,  w,  d,  theta 
DOUBLE  PRECISION  alpha,  h,  u 
DOUBLE  PRECISION  beta 

DOUBLE  PRECISION  ya,  yc,  phi,  a,  b,  c,  xO,  z,  xa,  xaprime,  zeta, 

&  zetaprime 

DOUBLE  PRECISION  fl,  f2,  f3,  f4,  fl  prime,  f2prime,  f3prime,  phia 
DOUBLE  PRECISION  phiaprime,  dphix,  dphixprime,  dphiz,  dphizprime 
DOUBLE  PRECISION  fn(l  1),  intz(21),  integralA,  integralC 
DOUBLE  PRECISION  mwbow,  mwstem,  mw 

**********  Constants 

pi=  4.D0  *DATAN(1  .DO) 

g=  32.2d0  !  ft/sA2 

rho=  62.4d0  /g  !  mass  density  of  water 

degrees2rads=  pi  /180.D0  !  conversion  factor:  degrees  to  radians 

**********  Barge  data 
WRITE(*,*)  'Enter  barge  wetted  length  (ft)' 

READ(*,*)L 

WRITE(*,*)  'Enter  barge  wetted  beam  (ft)’ 

READ(  *,*)  w 

WRITE(*,*)  'Enter  barge  draft  (ft)' 

READ(*,*)  d 

WRITE(*,*)  'Enter  bow/stem  rake  angle  (degrees)' 

READ(*,*)  theta 

theta=  theta  *degrees2rads  !  bow/stem  rake  angle,  rads 

**********  waterway  data 
WRITE(*,*)  'Enter  barge  speed  (ft/s)’ 

READ(  V)  u 

WRITE(*,*)  'Enter  water  depth  (ft)' 

READ(  V)  h 

WRITE(*,*)  'Enter  distance  from  bow  comer  to  wall  (ft)' 

READ(V)  ya 

WRITE(*,*)  'Enter  angle  of  approach  to  wall  (degrees)' 

READ(*,*)  alpha 

alpha=  alpha  *degrees2rads  !  wall  approach  angle,  rads 

**********  parametric  coefficients  in  the  analysis 

are  as  follows: 
yc=  ya  +L  *DSIN(alpha) 

phi=  0.5673d0  *(pi  +LdO  /DTAN(theta) )  *u**3  /g 
a=  u/  phi 
c=  pi  /h 


A2 
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b=  DSQRT(  DABS(-a*a  +c*c)  ) 
beta=  -b  *w  /2.d0 
x0=  pi  *u**2  /(2.d0*g) 

C  NOTE:  In  the  following,  prime  refers  to  the  image  barge 

C  **********  Bow  calculations: 
zeta=  -w/2 

C  Vertical  iteration  loop 
DO  10  j=  1,21 
z =  -d 

C  Transverse  iteration  loop 
DO  20  k=  1,11 
xa=  xO  +z  /DTAN(theta) 
xaprime=  xa  *(  2.d0  *DCOS(alpha)  -l.dO)  + 

&  w  *DSIN(alpha)  *DCOS(alpha)  +2.d0  *ya  *DSIN  (alpha) 
zetaprime=  zeta  *(  2.d0  *  DC  OS  (alpha)  -  l.dO)  - 
&  xa  *DSIN(2.D0*alpha)  +w  *(DCOS(alpha))**2  + 

&  2. dO  *ya  *DCOS(alpha) 

NOTE:  In  the  below  functions,  the  trigonometric  substitution 
l/COSH(x)  =  SECH(x)  =  SQRT(1  -(TANH(x))**2) 
is  used  (which  follows  from  the  identity 

(SECH(x))**2  +(TANH(x))**2  =  1 
This  is  done  because  COSH(x)  blows  up  as  x  becomes  large. 
fl=  DSQRT(  l.dO  -(DTANH(a  *xa))**2  ) 
f2=  DTANH(a  *xa) 
f3=  l.dO 
f4=  DCOS(c*z) 

flprime=  DSQRT(  l.dO  -(DTANH(a  *xaprime))**2  ) 
f2prime=  DTANH(a  *xaprime) 

f3prime=  DSQRT(  l.dO  -(DTANH(DABS(b  *zetaprime)  +beta))**2  ) 
phia=  -phi  *fl  *f3  *f4 
phiaprime=  -phi  *fl  prime  *f3  prime  *f4 
dphix=  a  *phi  *fl  *f2  *f3  *f4 
dphixprime=  a  *phi  *fl prime  *f2prime  *f3prime  *f4 
dphiz=  c  *phi  *fl  *f3  *DSIN(c*z) 
dphizprime=  c  *phi  *fl  prime  *f3 prime  *DSIN(c*z) 
fn(k)=  (phia  +phiaprime)  *(  (dphix  +dphixprime)  * 

&  DSIN(theta)  -(dphiz  +dphizprime)  * 

&  DCOS(theta)/DTAN(theta) ) 
z=  z+d  /lO.dO 
20  CONTINUE 

intz(j)=  (d/lO.dO)  *(  fn(l)  +4.d0*fn(2)  +2.d0*fn(3)  + 

&  4.d0*fn(4)  + 

&  2.d0*fn(5)  +4.d0*fn(6)  +2.d0*fn(7)  +4.d0*fn(8)  +2.d0*fn(9)  + 

&  4.d0*fn(10)  +fn(l  1)  ) 
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10  CONTINUE 


integralA=  (w/20.d0)  *(  intz(l)  +4.*intz(2)  +2.d0*intz(3)  + 

&  4.d0*intz(4)  +2.d0*intz(5)  +4.d0*intz(6)  +2.d0*intz(7)  + 

&  4.d0*intz(8)  +2.d0*intz(9)  + 

&  4.d0*intz(  1 0)  +2.d0*intz(  1 1 )  +4.d0*intz(  1 2)  +2.d0*  intz(  13)  + 

&  4.d0*intz(  14)  +2.d0*intz(  1 5)  +4.d0*intz(  1 6)  +2.d0*intz(  1 7)  + 

&  4.d0*intz(18)  +2.d0*intz(19)  +4.d0*intz(20)  +intz(21)  ) 

mwbow=  DABS(  rho  *  integral  A  /u**2  ) 

C  **********  Stem  calculations: 

zeta=  -w  /2.d0 

C  Vertical  iteration  loop 

DO  30  j=  1,21 
z=  -d 

C  Transverse  iteration  loop 
DO  40  k=  1,11 
xa=  xO  +z  /DTAN(theta) 
xaprime=  xa  *(  2.d0  *DCOS(alpha) -l.dO)  + 

&  w  *DSIN(alpha)  *DCOS(alpha)  +2.d0  *yc  *DSIN(alpha) 
zetaprime=  zeta  *(  2.D0  *DCOS(alpha)  -  l.dO)  - 
&  xa  *DSIN(2.D0*alpha)  +w  *(DCOS(alpha))**2  + 

&  2.d0  *yc  *DCOS(alpha) 

NOTE:  In  the  below  functions,  the  trigonometric  substitution 
l/COSH(x)  =  SECH(x)  =  SQRT(1  -(TANH(x))**2) 
is  used  (which  follows  from  the  identity 
(SECH(x))**2  +(TANH(x))**2  =  1 
This  is  done  because  COSH(x)  blows  up  as  x  becomes  large, 
fl  =  DSQRT(  1  .d0  -(DTANH(a  *xa))*  *2  ) 
f2=  DTANH(a  +xa) 
f3=  l.dO 
f4=  DCOS(c*z) 

flprime=  DSQRT(  l.dO  -(DTANH(a  *xaprime))**2  ) 
f2prime=  DTANH(a  *xaprime) 

f3prime=  DSQRT(  1  ,d0  -(DTANH(  DABS(b  *zetaprime)  +beta))**2  ) 
phia=  -phi  *fl  *f3  *f4 
phiaprime=  -phi  *fl  prime  *f3 prime  *f4 
dphix=  a  *phi  *fl  *f2  *f3  *f4 
dphixprime=  a  *phi  *fl  prime  *f2prime  *f3prime  *f4 
dphiz=  c  *phi  *fl  *f3  *DSIN(c*z) 
dphizprime=  c  *phi  *fl  prime  *f3prime  *DSIN(c*z) 
fn(k)=  (phia  +phiaprime)  *(  (dphix  +dphixprime)  * 

&  DSIN(theta)  -(dphiz  +dphizprime)  * 

&  DCOS(theta)/DTAN(theta) ) 
z =  z+d  /lO.dO 
40  CONTINUE 
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intz(j)=  (d/lO.dO)  *(  fn(l)  +4.d0*fo(2)  +2.d0*fn(3)  +4.d0*fn(4)  + 
&  2.d0*fh(5)  +4.d0*fn(6)  +2.d0*fn(7)  +4.d0*fn(8)  +2.d0*fn(9)  + 

&  4.d0*fn(10)  +fn(l  1) ) 

zeta=  zeta  +w/20.d0 
30  CONTINUE 

integralC=  (w/20.d0)  *(  intz(l)  +4.d0*intz(2)  +2.d0*intz(3)  + 

&  4.d0*intz(4)  + 

&  2.d0*intz(5)  +4.d0*intz(6)  +2.d0*intz(7)  +4.d0*intz(8)  + 

&  2.d0*intz(9)  + 

&  4.d0*intz(10)  +2.d0*intz(l  1)  +4.d0*intz(12)  +2.d0*intz(13)  + 

&  4.d0*intz(14)  +2.d0*intz(15)  +4.d0*intz(16)  +2.d0*intz(17)  + 

&  4.d0*intz(18)  +2.d0*intz(19)  +4.d0*intz(20)  +intz(21)  ) 

mwstem=  DABS(  rho  *integralC  /u**2  ) 

mw=  mwbow  +mwstem 

WRITE(*,*)  'Added  weight  =',  raw  *g  /2000.d0,  ’  tons.' 

STOP 

END 
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This  appendix  presents  sample  calculations  using  the  computer  code  presented 
in  Appendix  A.  Two  sample  cases  are  input,  and  the  corresponding  output  is 
listed. 

The  first  case  is  for  the  loaded  JOHB  abutting  the  wall  at  a  right  angle: 

Enter  barge  wetted  length  (ft) 

195. 

Enter  barge  wetted  beam  (ft) 

35. 

Enter  barge  draft  (ft) 

10. 

Enter  bow/stern  rake  angle  (degrees) 

45. 

Enter  barge  speed  (ft/s) 

3. 

Enter  water  depth  (ft) 

18. 

Enter  distance  from  bow  corner  to  wall  (ft) 

0. 

Enter  angle  of  approach  to  wall  (degrees) 

90. 

Added  weight  =  3.07247498044177  tons. 

The  second  case  is  for  the  loaded  JOHB  traveling  parallel  to  the  wall  in  open 
waters: 

Enter  barge  wetted  length  (ft) 

195. 

Enter  barge  wetted  beam  (ft) 

35. 

Enter  barge  draft  (ft) 

10. 
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Enter  bow/stern  rake  angle  (degrees) 

45. 

Enter  barge  speed  (ft/s) 

3. 

Enter  water  depth  (ft) 

18. 

Enter  distance  from  bow  corner  to  wall  (ft) 

70. 

Enter  angle  of  approach  to  wall  (degrees) 

0. 

Added  weight  =  3.14077173417141  tons. 
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